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Summary 

A self-consistent set of values have been computed for enthalpy, total specific heat at constant 
pressure, compressibility factor, viscosity, total thermal conductivity, and total Prandtl number 
of equilibrium air from 500 to 30 000 K over a pressure range of 10 4 to 10 2 atm. Previous results 
for most of these properties are limited to a temperature of 15000 K. The mixture values are 
calculated from the transport and thermodynamic properties of the individual species provided 
in a recent study by the present authors. The concentrations of the individual species, required 
in the mixture relations, are obtained from a free-energy-minimization calculation procedure. 
Calculations are based on an 11-species air model (O 2 , N 2 , O, N, NO, 0 + , N + , NO + , 0 ++ , 
N++ j e~ ). For pressures less than 10 — 2 atm and temperatures of about 15 000 K and 
greater, the concentrations of N++ and 0 ++ become important, and consequently they are 
included in the calculations for determining the various properties. The differences in the 
computed values given herein and those obtained by other researchers (which are based on 
temperatures to 15 000 K) may be attributed to the differences in the properties employed for 
the individual species, mixing laws used for the mixture, and (to a lesser degree) the number 
of chemical species considered in the mixture. Finally, all the computed properties are curve 
fit as a function of temperature at a constant value of the pressure. These curve fits reproduce 
the computed values to within 5 percent for the entire temperature range considered at specific 
pressures and provide an efficient means for computing the flow-field properties of equilibrium 
air, provided the elemental composition remains constant at 0.24 for oxygen and 0.76 for nitrogen 
by mass. 

Introduction 

Equilibrium air chemistry calculations can be performed either by using the free-energy- 
minimization procedure (ref. 1) or by using the equilibrium constant method (refs. 2 and 
3). Flow-field calculations based on either of these approaches (refs. 4 and 5) can be quite 
time consuming; however, the assumption of chemical equilibrium without elemental separation 
does allow the governing equations to be written in a form such that the individual species 
concentrations are not required explicitly. In this format, the specific heat at constant pressure 
Cp , the thermal conductivity A, and the Prandtl number iVp r are defined as total properties 
(refs. 6, 7, and 8). Hansen (ref. 7) has provided these properties in a tabular form for 
temperatures up to 15 000 K and pressures from 10 4 to 10 2 atm using a seven-species air 
model (N 2 , 0 2 , N, O, N+, 0+, and e“). Peng and Pindroh (ref. 9) provide tabulated values 
of total thermal conductivity for a nine-species air model (N 2 , O 2 , N, O, NO, N + , 0 + , NO + , 
and e~); however, no values are given for the total specific heat or the Prandtl number. The 
tabulated values of reference 9 are for temperatures up to 15000 K and pressures of 10 -2 to 
10 2 atm. Both reference 7 and reference 9 have incorporated a lower order electron pressure 
correction (compared with the one used in ref. 8) in the transport properties of ionic species to 
obtain the mixture values for equilibrium air. The thermodynamic properties provided in these 
two references and in reference 10 1 are based on an equilibrium air model that did not include 
N ++ and 0 ++ among its constituent species. These two species are not significant at pressures 
of 10” 4 to 10 2 atm for temperatures less than 15 000 K. However, for higher temperatures (such 
as those considered in ref. 10), these two species (especially N ++ ) become quite significant. (See 
fig. 1.) 

Since curve fits to the various properties permit efficient computation for flow-field analyses 
(ref. 11), it is desirable to have the curve fits in preference to the tabulated values. Further, 
for accurate calculations, it is imperative that these curve fits be obtained to values that are 


1 Vinokur and his colleagues have recently included doubly ionized atomic species in their calculations of the thermody- 
namic properties. These results, however, have not been published. 


computed in a self-consistent manner from the same set of data and other physical properties. 
Srinivasan et al. (refs. 12 and 13) recently obtained curve fits to the tabulated transport 
properties of reference 9 and the thermodynamic properties of reference 14. No curve fits are 
readily available for tabulated values of Hansen (ref. 7). 

This report presents accurate curve fits for the computed values of enthalpy, total specific 
heat, compressibility factor, viscosity, and total values of thermal conductivity and Prandtl 
number of equilibrium air over a more extensive range of parameters than had been available. 
This has been accomplished by using both the thermodynamic and the transport properties 
of the individual species (in the form of curve fits) from reference 8 along with a free-energy- 
minimization procedure (ref. 1) for computing the equilibrium composition. Figure 1 shows the 
calculated composition of equilibrium air (used in the calculations of the properties obtained 
herein) as a function of temperature for pressures of 10 -4 , 10^, and 10 2 atm. As shown in 
figure 1(b), the doubly ionized atomic species are dominant at a pressure of 10“ 4 atm for 
temperatures of about 15 000 K and greater. Their appearance and dominance are delayed to 
higher temperatures at higher pressures. (See figs. 1(d) and 1(f).) With the individual species 
composition and properties known, the mixture properties for an 11-species air model (N 2 , O 2 , 
NO, N, O, N 4 ", 0 + , NO 4 ", N^ 4 ", O 4-4 *, and e“) are then calculated. These computed mixture 
properties are curve fit for temperatures of 500 to 30000 K and for specific pressures of 10~ 4 
to 10 2 atm. Calculations are limited to 25000 K at 10 -4 atm and 28 000 K at 10~ 3 atm to 
avoid the necessity of including contributions from triply ionized species. Tabulated values of 
these curve-fit coefficients are given herein and computer subroutines to evaluate the various 
properties using these coefficients are provided in reference 15. These values and codes may be 
obtained for a fee from the following: 

COSMIC 

Computer Services Annex 

University of Georgia 

Athens, GA 30602 
(404) 542-3265 

Request the code by the designation LAR 14760. This code is written in FORTRAN 77 for use 
on computers with a FORTRAN compiler. 

Detailed comparisons of the curve-fit data with computed results are presented to demonstrate 
the accuracy of the obtained curve fits. Further, the presently computed values are compared 
with existing data in the literature. Comparisons are also included to show the effect of the 
number of species included in an air model on the equilibrium properties. 

The curve fits provided herein can only be used for equilibrium air calculations without 
ablation or with or without air injection, that is, with constant elemental composition through 
the flow field. Values of 0.24 and 0.76 for the elemental composition by mass of oxygen and 
nitrogen, respectively, are employed in this work. 

The properties are obtained herein under the assumption of charge neutrality for the ionized 
air. Therefore, the associated electric field effects are neglected. This assumption can generally 
be made for an ionized gas mixture in thermochemical equilibrium. Further, the transport 
properties are computed from the approximate formulas of references 8, 9, 16, 17, and 18 
instead of the complete kinetic theory expressions of reference 19. These approximate formulas, 
although still somewhat complicated, are much simpler to use than the equations given by the 
first Chapman-Enskog approximation (ref. 19) and give results within a few percent of those of 
reference 19 for unionized air through temperatures of 8000 K (ref. 20). Since the accuracies 

of collision integrals fi d ' and fi-d ' at high temperatures are known only to about 10 to 


2 


20 percent at best, the error introduced in our calculations by the use of approximate formulas in 
place of the complete Chapman- Enskog approximation may be, therefore, relatively insignificant. 
For more exact treatment, however, the formulas given in reference 8 (based on the relations 
developed by Yos, ref. 21) may be used to compute the mixture properties. 

The accuracy of the basic calculations carried out herein for the transport properties of 
partially ionized air is difficult to evaluate without reliable experimental data for the temperature 
and pressure ranges considered in this report. Detailed comparisons, however, are provided with 
the existing calculations of Hansen (ref. 7) and of Peng and Pindroh (ref. 9). Devoto and Bauder 
(ref. 22) have obtained data for the thermal conductivity of equilibrium air from electric arc 
measurements at a pressure of 1 atm and temperatures of 8000 to 13000 K. Their data suggest 
somewhat higher values for the thermal conductivity when appreciable ionization is present. The 
electric arc data, however, are not considered to be very accurate because the thermodynamic 
state of the gas may not be well defined, with some of the internal degrees of freedom partially 
frozen. This lack of accuracy is reflected in the large scatter of Devoto’s data. 

Finally, the effects of Coulomb collisions are not obtained very accurately with the first 
Chapman-Enskog approximation. In the limiting case of a fully ionized gas, the transport 
properties calculated with this approximation differ from those obtained with an accurate 
solution of the Boltzmann equation (refs. 23 and 24) by about a factor of 2. References 23 
and 24 have also suggested approximations that account for effects of Coulomb collisions more 
accurately. But a complete solution of the problem is not given for a gas with an arbitrary degree 
of ionization. In the present calculations, an interpolation procedure is used for the region 
of partial ionization and involves the approximate formulas (used in the transport properties 
calculations and given later) using effective Coulomb collision cross sections. These cross sections 
are chosen (similar to ref. 17) to make the calculated transport properties agree as closely as 
possible with the results of Spitzer and Harm (ref. 23) for the fully ionized case. 


Symbols and Abbreviations 


A C P ’ B c p > c c p - D c v > E c p 

A h' B h’ C / ii B h' E h 

A Ki b k, C K> d k> e k 

A V r , £p r , Cp r , Dp r , E P r> F Pr 

A Z> B Zi C Zi D z > E z 
A fx, Bp, Cp, Dp, Ep , F 'n 
Ci 

c p 

c Pf 


curve- fit coefficients for specific heat at constant 
pressure Cp (eqs. (14a) and (14b)) 

curve-fit coefficients for enthalpy h (eqs. (13a) and 
(13b)) 

curve-fit coefficients for total thermal conductivity K 
(eqs. (30a) and (30b)) 

curve-fit coefficients for Prandtl number Np r 
(eq. (31a)) 

curve-fit coefficients for compressibility Z (eq. (19a)) 
curve-fit coefficients for viscosity fi (eq. (22a)) 
mass fraction of species z, pj p 
specific heat at constant pressure, 

C p = Cpj T Cp r — (^37^ ’ cal/g-K 

frozen specific heat at constant pressure, 

NS 

£ CiC pii , cal/ g-K 
i— 1 
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Cp r 

reactive specific heat at constant pressure, 
|, j (^) p , ca ,/ g . K 

Cp,i 

specific heat at constant pressure of species z, 

(w ); -‘/ g -K 

Dij 

binary diffusion coefficient, cm 2 /sec 

Dij 

multicomponent diffusion coefficient, cm 2 /sec 

Dim. 

effective binary diffusion coefficient (eq. (4b)), cm 2 /sec 

Dim 

effective binary diffusion coefficient (eq. (5b)), cm 2 /sec 

H 

u 2 

total enthalpy of mixture, h + — , kcal/g 

h 

NS 

enthalpy of mixture, ^ Cjh^ kcal/g 

z= 1 

hi 

enthalpy of species z, kcal/g 

Ji 

diffusion mass flux of species z, g/cm 2 -sec 

3 

zero for plane flow and one for axisymmetric flow 

K 

total thermal conductivity, K j 4- K r , cal/cm-sec-K 

K f 

frozen thermal conductivity, cal/cm-sec-K 

Kmt 

internal component of frozen thermal conductivity, 
cal/cm-sec-K 

K r 

reaction component of total thermal conductivity, 
cal/cm-sec-K 

K s 

thermal conductivity from Sutherland’s law (eq. (33)), 
cal/cm-sec-K 

K tr 

translational component of frozen thermal conductiv- 
ity, cal/cm-sec-K 

k 

Boltzmann’s constant, 1.38066 x 10 - ^ erg/K 

M 

molecular weight of mixture, g/g-mole 

Mi 

molecular weight of species z, g/g-mole 

M 0 

molecular weight of undissociated air, 28.96 g/g-mole 

n a 

Avogadro’s number, 6.0221 x 10 23 molecules/g-mole 

iVpr 

total Prandtl number, Cpfi/K 

NPrJ 

frozen Prandtl number, Cpj(i/Kj 

n 

coordinate measured normal to body, cm 

V 

pressure, atm 

Q 

wall heat transfer rate, cal/cm 2 -sec 
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-^univ 

T 


S 

T 

u 

V 


Ul; 


Z 




e 


K 

V 

Vs 

P 

Pi 

Po 

n|5 ,1} 


n< 2 ' 2 > 


Subscripts: 

e 

i 

3 

l 

m 

V 

Abbreviations: 

NIR 

NS 


universal gas constant, 1.987 cal/g-mole-K 

radius measured from axis of symmetry to a point on 
body surface, cm 

coordinate measured along body, cm 
temperature, K 

velocity component tangent to body surface, cm/sec 
velocity component normal to body surface, cm/sec 
mass rate of formation of species i, g/cm^-sec 
mole fraction of species i 
compressibility factor 

defined by equation (26), cm-sec 

defined by equation (21), cm-sec 
local body angle measured from body axis 
body curvature, cm - * 
viscosity, poise 

viscosity from Sutherland’s law (eq. (32)), poise 

density, g/cm’* 

density of species g/cnV* 

density of undissociated air, g/cm^ 

average collision cross section (used for thermal 
conductivity) for collisions between species i and j, 
A 2 (lA — 10 - ^ cm) 

average collision cross section (used for viscosity and 
thermal conductivity) for collisions between species i 
and j, A 2 

local chemical equilibrium 

species i 

species j 

species l 

mixture 

constant pressure 

total number of independent reactions 
total number of chemical species 
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Background 

This section addresses the usefulness of employing the total properties in place of the frozen 
values to solve the flow-field equations under chemical-equilibrium conditions. The energy 
equation for a chemically reacting multicomponent gas mixture can be written in terms of 
frozen values of thermal conductivity and specific heat at constant pressure. (See ref. 4, for 
example.) However, if the system remains in equilibrium and the elemental composition is 
constant throughout the flow (i.e., there is no foreign gas injection or elemental separation 2 ), 
then the equation can be written in terms of total thermal conductivity and total specific heat 
at constant pressure. (See appendix A.) In such a formulation, the species production term 
ihj, the species diffusion mass flux J ; , and species concentration C t do not appear explicitly in 
the energy equation but instead are contained in the definitions of total specific heat and total 
thermal conductivity. 

Now the total specific heat at constant pressure C p is defined as follows (see appendix A and 
refs. 7, 8, and 25): 


r ( dh \ 

P \dTj p 

(la) 

NS NS \ 

= £ c ; c p,; + £M a#) 
i= 1 i=l V / P 

(lb) 

— C P j + C Pr 

(lc) 


The first term of equations (lb) and (lc) is the frozen component of the specific heat and the 
second term represents the contribution from chemical reactions. Similarly, the total thermal 
conductivity K contains the frozen and reactive components and may be defined as (see refs. 6, 
7, 8, 17, and 25) 

K = K f + K r (2a) 


or 


K = K f - (^hiJ^idT/dn) 


(2b) 


where Kj consists of the translational K iv and internal AT int components, which are defined 
subsequently (eqs. (24a) and (25)). 

The diffusion mass flux of species J j, in equation (2b) is given by (refs. 19, 26, and 27) 


NS 


Ji — 2 5Z A fiMjDjj 


M P i 


dx j 

dn 


( 3 ) 


for multicomponent diffusion with pressure and thermal diffusion neglected. For a mixture of 
more than two species, the multicomponent diffusion coefficients D^j depend not only on the 
properties of the individual species i and j but also on the mixture composition and on the 
properties of all species in the mixture (ref. 28). These multicomponent coefficients can be 
calculated with the expressions given in reference 19 from the mixture composition and the 
ordinary binary diffusion coefficients Djj given in reference 8. A simpler form of equation (3) 

2 Elemental separation under thermochemical equilibrium conditions is generally negligible unless there is a foreign gas 
injection or ablation at the surface of a body. 
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can be obtained with the assumption that all the species in the mixture move with the same 
velocity or are stationary. With this assumption, equation (3) in terms of the mass-fraction 
gradients becomes (refs. 9, 25, and 26) 


with 


Jj = -pD im ^ 
1 H 1171 dn 

1 — Xi 




53 X j/D t j 

or, in terms of the mole- fraction gradients (refs. 29 and 30), 

T - - Ti 
Ji ~ pUim dn 


(4a) 

(4b) 


(5a) 


with 


^im — 


{Mj/ M) [l - (A/^/M) Xj\ 


5^ X j/ D ij 
3 ¥* 


(5b) 


Use of any of the above formulations to obtain total thermal conductivity K requires the 
solution of equation (2b) simultaneously with the flow-field governing equations for a particular 
flow problem. Under the chemical equilibrium condition, however, the chemical composition 
almost comes into equilibrium with the local temperature (ref. 6). For such a case, therefore, 
one can write 


dC t \ dCi/dn 
dT ) e ~ dT/dn 


(6) 


{ dxi \ dxi/dn 

\df) e = WJ&n 


where the subscript e implies local chemical equilibrium. Using equations (6) and (7) along with 
equations (3), (4a), and (5a) in equation (2b), one obtains 


NS NS /F> N 

= %+ 3 EEW¥ij # 

M i= lj=l V°N/e 


(8a) 



K = K f + P '£D im h l (|£) 
*= 1 V J e 

(8b) 

NS s 

K^Kf + p^Dimhil^) 

i=l \ U1 / e 

(8c) 


Evaluation of K from either the exact equation (8a) or the approximate equation (8b) or (8c) 
does not require any details of the flow-field geometry. Obviously, equation (8b) or (8c) is 
considerably simpler to use than equation 8(a) for the evaluation of K. These equations imply 
that the total thermal conductivity is being computed at a constant pressure. 
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If the energy equation is formulated in terms of total enthalpy instead of temperature, then 
it can be shown that the equation can be written in a simplified form (see appendix A) with the 
use of one additional total quantity, the total Prandtl number iVp r , defined as 


iVp r = C p fi/K ■ (9) 

where the total values of C p and K are obtained from equations (1) and (2), respectively. 

Thus, with the introduction of total properties, it is possible to write the energy equation in 
terms of either temperature or total enthalpy and the resulting equations have the same form as 
the corresponding perfect-gas equations. Consequently, the equations for an equilibrium reacting 
gas mixture can be solved in a manner similar to that for a perfect gas, provided the curve fits 
to the total values of C p , K, and 7Vp r are available for the ranges of temperatures and pressures 
of interest. Such curve fits have been obtained in this study. As also mentioned previously, these 
curve fits can be used only for equilibrium air calculations where the elemental composition is 
constant throughout the flow field. 

No curve fits are given here for the frozen values of Cp, K, and lVp r . This is because the 
energy equation formulated in terms of these quantitites (see eqs. (A2) or (A4)) involves gradients 
of the species concentrations. This formulation, therefore, requires an explicit solution for the 
species concentrations in the flow field. Once the concentration (and temperature) distributions 
are available, the frozen values of Cp, K, and lVp r can easily be obtained from the individual 
species properties. Therefore, there is no need to obtain curve fits for the frozen values. 

Analysis 

Curve-Fit Temperature Intervals, Averaging, and Interpolation Procedures 

The order of the curve-fit polynomials and the number of temperature intervals for the curve 
fits have been selected to achieve a 5-percent overall accuracy of the curve-fit values compared 
with the computed values. Also, less computational effort is needed in computing the properties 
in flow-field calculations if fewer temperature intervals (or ranges) are used with a lower order 
polynomial. The entire temperature range of 500 to 30000 K is broken down to two or more 
intervals with these criteria. Further, to ensure a smooth variation of the curve-fit properties 
over the complete temperature range, values of the curve-fit coefficients are linearly averaged 
at the temperature interval boundaries over a specified overlap region. The overlap region is 
equally spread on the two sides at the temperature boundaries and is 500 K for temperatures 
less than 20 000 K and 1000 K for higher temperatures. The averaging procedure for a coefficient 
in the overlap region can be written as 


A <j> = (! ~ a ) A % + aA $ 


(10a) 


where 


a = 


( T-Tk 

[^~ t b 


(10b) 


Here, Tg and are the lower and upper temperature limits, respectively, of the overlap region. 
Coefficient A ( . } is one of the curve-fit coefficients for flow-field property <f>- The superscripts U and 
L on denote coefficients of the curve fits from the upper and lower sides of the temperature 

boundaries, respectively. With these definitions, the curve-fit coefficient A ^ is replaced by A ( p 
from equation (10a) when the temperature at which the property is to be evaluated is within an 
overlap region. An example subroutine, which evaluates the curve fits and performs the linear 
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averaging, is presented in appendix B. This subroutine can be easily modified to suit the user’s 
requirements. 

Since the curve fits are obtained for a range of constant pressures (10 4 to 10^ atm) in 
steps of a factor of 10, an interpolation procedure is required to obtain values of the flow-field 
properties at intermediate pressures. A simple logarithmic interpolation procedure is suggested 
here for this purpose. With this procedure, a property </>* may be interpolated at pressure p* 
(for a specified temperature) from the following relation: 

log^* = _ iogp| ) (log P* - log Pi) + lo g^l ( U ) 

where <\> 2 and (j> \ are the values of property 0 evaluated for a specified temperature T at pressures 
P 2 and p \ , respectively, such that p\ < p* < P2- Other details of the interpolation procedure are 
given in reference 15. This reference also provides an alternate interpolation procedure, which 
gives smaller interpolation error as compared to equation (11). 

The curve fits provided herein are in terms of temperature. These may be easily converted 
to be in terms of enthalpy by using the relation between temperature and enthalpy provided by 
the curve fits in the next section. 

Details of Property Computations and Curve Fits 

Enthalpy . The specific enthalpy for equilibrium air can be computed from 

NS 

h = J2 c i h i ( 12 ) 

z=l 

where h ^ is the specific enthalpy of individual species obtained from the curve fits of reference 8 
and C l is the concentration of individual species obtained from a free-energy-minimization 
calculation procedure (ref. 1). These computed values have been curve fit by employing the 
following curve-fit expression to obtain enthalpy as a function of temperature at constant 


pressure: 

h = exp (A/jX 4 + B hX 3 + c hX 2 + D hX + E h ) 

(13a) 

or 

In h = A/ t x 4 + B hX 3 + E hX 2 + B hX + E h 

(13b) 

where 

X = ln(T / 1 0 000) 

(13c) 


The coefficients A^, Bf v C / l , D^, and have been evaluated from the computed values 
with a least-squares curve-fit technique and are given in table I. These coefficients are given 
for several temperature ranges at a specified pressure. As mentioned previously, the number of 
temperature ranges is established to obtain a smooth and accurate curve fit to the computed 
values, and generally a larger number of ranges are needed at lower pressures. 

Figure 2 shows the curve-fit and computed values for enthalpy at selected pressure values. 
All the computed values are not shown in this and other figures, especially at temperatures of 
less than 2000 K. Values from the curve-fit equation (13a) are within 4 percent of the values 
computed from equation (12) for the seven specified pressures in the range 10 -4 to 10 2 atm and 
the entire temperature range (500 to 30 000 K) considered herein. 
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Total specific heat at constant pressure. The total specific heat at constant pressure can 
be computed from equation (lb). The species concentration Q appearing in this equation is 
obtained from the free-energy-minimization calculation method of reference 1, and (dCJdT) p is 
evaluated numerically by differentiating the data from this calculation. The contribution of the 
chemical reactions to the specific heat (involving {OC,J dT) p ) becomes neligible for temperatures 
less than 2000 K at pressures around 10® atm. 

The following expression is employed to curve fit the values of C p computed from equa- 
tion (lb): 

C p = exp (Ac p X A + B Cp x 3 + C Cp X 2 + D C p X + E C P ) ( 14a ) 

or 

In C p = A Cp x 4 + B Cp x 3 + Cc p X 2 + D C p X + E Cp (14b) 

with 

X = ln(T/10000) (14c) 

The tabulated values of coefficients Aq^ Cq^ and Eq are provided in table II. 

Figure 3 shows the computed and curve-fit values. The curve-fit values are within 4 percent of 
the computed values for the entire temperature range and pressure values considered herein. 


Compressibility factor . The density of equilibrium air can be obtained from the equation 
of state: 

p = (pM/R unW T) x 2.4218 x 10 -2 (15) 

where the factor 2.4218 x 10 — 2 appears because of the units of thermodynamic variables and 
the universal gas constant employed here (see list of symbols for the units) and the mixture 
molecular weight M is given by 


M 


{Ci/M,) 


The density of undissociated air p 0 can be computed from 


(16) 


p 0 = (pA/o/i? un j v T) 2.4218 x 10 -2 


(17) 


The compressibility factor Z, which is equal to the ratio of the molecular weight of undissociated 
air to the mean molecular weight, can be obtained from equations (15) and (17) through the 
relation 


Mq _ Po 

W ” ~p 


(18) 


The computed values of the compressibility factor Z have been curve fit by employing the 
following expression at constant pressure: 


Z = A z + B z x + C z x 2 + D z x 3 + E z x 1 


(19a) 


where 

X = T/1000 (19b) 

The coefficients appearing in equation (19a) are given in table III. Accuracy of the curve fit of 
equation (19a) compared with the computed values from equations (15), (17), and (18) is shown 
in figure 4. The curve-fit values are within 3 percent of the computed values for the pressure 
range of 10 -4 to 10 2 atm. 
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Viscosity. Viscosity has been computed from the following expression (see refs. 8 and 17): 


NS 


*■ = £ 


(Mi/N A )x i 


: — ' y^ NS x . a (2) 

*=1 


( 20 ) 


(2) 

where A-- is defined as 




A^ = —(1,5460 x 10 -20 ) 
l 3 5 


2M l M J 


1/2 


^ ^univT (Vi + 


*n<f> 


( 21 ) 


In the above equation, is the average collision cross section for viscosity and is provided 

in curve-fit form a s a function of temperature in reference 8 for the various species interaction 
pairs. 

The following expression has been used to curve fit the viscosity values computed from 
equation (20): 

M = Ap + B^x + + £>/iX 3 + £ M X 4 + ^/iX 5 (22a) 

where 

X = T/1000 (22b) 

The polynomial coefficients of equation (22a) are given in table IV. The computed and curve-fit 
values of the viscosity are shown in figure 5 for different pressures. In general, the curve-fit. 
values are within 4 percent of the values computed from equation (20). The accuracy improves 
with increasing pressure, similar to the other thermodynamic and transport properties. 


Total thermal conductivity . The total thermal conductivity K for equilibrium air is 
obtained from equations (2). The frozen part in these equations is obtained from the translational 
and internal components: 

K f = K tr + K mt (23) 

The translational contribution in the above equation is computed from the relation (see refs. 8 
and 17) 


K tx = (2.3901 x 10- 8 )^£ 


NS 


4 y* NS x . A (2) 
i= 1 l j=\ a ij x 3 i-i 


(24a) 


where a| 2 ^ is given by equation (21) and a t j is defined as 


a ij — 1 + 


1 - 1 

[Mi/MjJ 

0.45 - 2.5^ 

1 (Mi/Mj)' 


1 + 


2 


(24b) 


The component of thermal conductivity resulting from the internal excitation energy of the 
molecules K mt is calculated from (refs. 8 and 17) 


NS 

K mt = (2.3901 x 10- 8 )fc£ 

i= 1 


'(gp.t/fiuniv) - (5/2)' 


y-NS 


x-A (1) 
l X 3 ij 


(25) 
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where A 


(i) 

ij 


is defined as 


1/2 

(26) 

and is the collision cross section for diffusion, which is provided in curve-fit form in 

reference 8. 

The reaction part of the thermal conductivity K r is defined in equation (8a) or, through 
some approximations, in equations (8b) and (8c). As shown for a single reaction in references 31 
and 32, the reaction part of equation (8a) is equivalent to 

K r = k X 10" 4 V univ r ) 2 ( 2 7) 

ti E|S E/Ji [vi ~ Vi) A{j> 

given in references 7, 8, 9, and 17. Compared with the definition of K r in equation (8a), 
equation (27) has been simplified with the expressions from reference 19 to eliminate the 
multicomponent diffusion coefficient Djj from the problem. Equation (27), therefore, can be 
evaluated with a greater degree of accuracy and has been used to evaluate K r . Equation (27) is 
equivalent to the reaction part of equation (8a) for those temperature and pressure conditions 
for which only a single chemical reaction is significant; however, there will be some error in 
this approximation for conditions in which more than one reaction is proceeding simultaneously. 
This error has been discussed, for example, in reference 20. 

In equation (27), NIR is the total number of independent reactions in the system and a[ ? 
represents the stoichiometric coefficients of components A j in the Ith chemical reaction written 
in the balanced form: 

NS 

= ° (28) 
i=l 

and Ah i is the heat of reaction per gram-mole for the Ith reaction given by 

NS 

Ah[ = y ^ xio ) (29) 

i~ 1 

The total thermal conductivity K , computed from equations (2a), (23), (24a), (25), and (27), 
is curve fit with the expression 

K = exp (. A k x 4 + B kX 3 + C K x 2 + &KX + Er) (30a) 

or 

In K = A k x 4 + Bkx 3 + CrX 2 + B>KX + E K (30b) 

where 

X = ln(T/10000) (30c) 

The coefficients appearing in equations (30) are given in table V. Figure 6 shows the computed 
and curve-fit values. The curve-fit values are within 4 percent of the computed values for the 
entire temperature range and pressure values considered herein. 


= |(1.546 x 10“ 20 ) 


2M t Mj 


L^univ^ Mi + Mj\ 
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Total Prandtl number. The total Prandtl number for equilibrium air has been computed 
from equation (9), with C p and p obtained from equations (1) and (20), respectively, and K 
calculated from equations (2a), (23), (24a), (25), and (27). The computed values of the Prandtl 
number have been curve fit with the relation 

iVpr = Ap v + Bp v x + CprX 2 + ^PrX 3 + E PrX 4 + ^Pr* 5 (31a) 

with 

X = T / 1000 (31b) 

The polynomial coefficients of equation (31a) are provided in table VI. Figure 7 shows the 
curve-fit values of equation (33a) compared with the computed values. The curve-fit values are 
within 4 percent of the computed values for the pressure and temperature ranges considered 
herein. 

Variation of the Properties With Pressure 

The curve- fit values for enthalpy, total specific heat at constant pressure, compressibility 
factor, viscosity, total thermal conductivity, and total Prandtl number are shown in figure 8. 
This figure gives the variation of these properties with pressure. For a given pressure, the 
temperature dependence of the various properties can be divided roughly into four portions 
corresponding to the reactions taking place. From 500 K to the temperature at which the 
curves begin to diverge, the mixture stays essentially inert. Divergence can be attributed to the 
dissociation of oxygen. The large change in properties at the next temperature portion occurs 
basically because of nitrogen dissociation. At the third temperature portion, single ionization 
of nitrogen and oxygen atoms causes the change in properties. Finally, double ionization of 
nitrogen and oxygen is responsible for the change in properties at the last temperature portion. 
These changes in the values of Cp, /i, K , and Vp r are reflected, generally, through maxima in 
their values. With the decrease in pressure these maxima (except for p) increase in sharpness 
as they shift to lower temperatures. 

Comparison With Existing Results 

A comparison is provided in this section of the presently computed results with those existing 
in the literature and commonly employed in various flow-field computations. 


Enthalpy. The computed dimensionless enthalpy values are compared with those of Hansen 
(ref. 7) in figure 9 for two pressures. As mentioned in the "‘Introduction,” Hansen’s values are 
available only through a temperature of 15 000 K for the pressure range of 10“^ to 10 2 atm. 
Agreement is good between the two results through this temperature range. The present results 
have been obtained with an 11-species air model (N 2 , 0 2 , N, O, NO, N+, 0 + , NO+, N++, 
0 ++ , and e ), whereas the results of reference 7 are based on a 7-species air model (N 2 , 0 2 , N, 
O, N+, 0 + , and e“). Obviously the small concentrations of NO+, N+ + , and 0+ + [C { < 10 -3 ) 
do not affect the equilibrium air enthalpy significantly through temperatures less than 15 000 K. 
For temperatures greater than 15 000 K and pressures less than 10~ 2 atm, however, N ++ and 
0 ++ concentrations become quite significant and the equilibrium air calculations need to include 
these two species, as shown in figures 10(a) and 10(b). A nine-species air model (N 2 , 0 2 , N, O, 
NO, N+, 0 + , NO + , and e ) gives enthalpy values similar to those obtained from a seven-species 
air model. The enthalpy values obtained from these two models differ considerably from those 
obtained from an 11-species air model at temperatures greater than 15 000 K. 


Total specific heat at constant pressure. The dimensionless values of the total specific 
heat at constant pressure are compared with those of Hansen (ref. 7) in figure 11. Agreement 


13 


is quite good except at low pressures and high temperatures, when the production of N ++ and 
0++ begins to dominate over other species. Clearly the specific heat calculations are affected 
more than the enthalpy calculations by the production of these two doubly charged ions. 


Compressibility factor. The compressibility factors computed presently and those given by 
Hansen (ref. 7) compare quite well for the entire pressure range (10~ 4 to 10 2 atm) through a 
temperature range of 15000 K, as shown in figure 12. Similar to the enthalpy values, differences 
in the number of chemical species used in the two calculations do not affect the compressibility 
factor to any noticeable degree for these pressure and temperature ranges. 


Viscosity. The presently computed viscosity values are compared with those of Hansen 
(ref. 7), Peng and Pindroh (ref. 9), and Svehla (ref. 33) in figure 13(a) for a pressure of 10° atm. 
Present results obtained with the 11-species air model compare well with the calculations of 
Peng and Pindroh (based on a 9-species air model) through a temperature of about 7500 K. 
The discrepancy between the two calculations at higher temperatures may be because of the 
differences in the mixing law and the cross sections used for charged-neutral and charged-charged 
particle interactions. Agreement of the present values with those of Hansen is quite poor for 
almost the entire temperature range shown in figure 13(a), except at temperatures less than 
about 1500 K. In fact, Hansen’s values are similar to those predicted from Sutherland’s law of 
viscosity 

T 3 / 2 

^» = 14584 x lo^fTTIius (32) 

for temperatures less than 5000 K. The agreement of these two values in the temperature range 
of 2000 to 5000 K is somewhat surprising in view of the fact that oxygen begins to dissociate 
at 2000 K and is fully dissociated at 5000 K. Sutherland’s law of viscosity is applicable only 
to undissociated air. All the computed values approach Sutherland’s values at temperatures 
less than 2000 K, as shown in figure 13(a). Also included in this figure are the viscosity values 
computed by Svehla (ref. 33) for unionized air. These values compare well with the present 
values and with those predicted by Peng and Pindroh through a temperature of about 4000 K. 

Figure 13(b) shows a comparison of the presently computed viscosity values with those of 
Hansen (ref. 7) and those from Sutherland’s law of viscosity (eq. (32)) at a pressure of 10~ 4 
atm. The disagreement with Hansen’s values is similar to that noted for a pressure of 10® atm 
(fig. 13(a)). Peng and Pindroh (ref. 9) have not provided viscosity values for p < 10 -2 atm and, 
therefore, their calculations are not included for comparison here. Similarly, Svehla (ref. 33) has 
provided computed values only for 10° atm and his values cannot be compared in figure 13(b). 

The differences in the various results shown herein, particularly at higher temperatures, may 
be attributed to the differences in the viscosities of the constituent air species (resulting from 
the differences in the species cross sections), mixing laws employed for the mixture calculations, 
and to a somewhat lesser extent, the number of chemical species considered in the air model. 


Total thermal conductivity. The present values of the total thermal conductivity are 
compared with those given by Hansen (ref. 7) and Peng and Pindroh (ref. 9) in figure 14(a) 
for a pressure of 10° atm. Figure 14(b) gives a comparison with Hansen’s (ref. 7) predictions 
for a pressure of 10~ 4 atm. As with viscosity, the total conductivity values of reference 9 
are not available for pressures less than 10 -2 atm. Further, the present thermal conductivity 
values agree better with the results of Peng and Pindroh than with those of Hansen through a 
temperature of about 9000 K, as shown in figure 14(a). With increasing ionization at higher 
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temperatures the two results begin to differ. Hansen’s results underpredict thermal conductivity 
through a temperature of about 9000 K and overpredict beyond 9500 K compared with the 
present results (as well as those of Peng and Pindroh). Also included in figures 14(a) and 14(b) 
are the values obtained from Sutherland’s law: 


r. T 3 / 2 

K s = 5.9776 x 1CT 6 - — — — 
T + 194.4 


(33) 


All the computed values approach the Sutherland values as the temperature decreases to 500 K, 
when only the rotational component of the internal energy of the molecular species is fully 
excited. For temperatures greater than 500 K, the vibrational component becomes excited and 
the dissociation of oxygen in air is initiated at approximately 2000 K. Therefore, Sutherland’s 
law (eq. (33)) should not be applicable for temperatures greater than about 500 K. 

The differences between the values for various computations shown in figures 14(a) and 14(b) 
exist for reasons similar to those given previously for the differences in viscosity values. 


Total Prandtl number. The calculated total Prandtl number values compare qualitatively 
with those of Hansen (ref. 7) for a temperature up to about 9000 K at a pressure of 10® atm, 
as shown in figure 15(a). However, the values differ quantitatively by as much as 22 percent. 
For temperatures greater than 9000 K, Hansen’s values are much lower than the present results, 
primarily because of the larger value of the thermal conductivity computed by Hansen. (See 
fig. 14(a).) Results from the two calculations compare much more favorably at a pressure of 
10“ 4 atm (as shown in fig. 15(b)) except for the temperature range of 5500 to 10500 K. These 
differences are related directly to the differences in the values of thermal conductivity shown in 
figure 14(b). 

Concluding Remarks 

A self-consistent set of values have been computed for enthalpy, total specific heat at constant 
pressure, compressibility factor, viscosity, total thermal conductivity, and total Prandtl number 
of equilibrium air from 500 to 30 000 K over a pressure range of 10“ 4 to 10 2 atm. Previous 
results for most of these properties have been limited to a temperature of 15 000 K. The mixture 
values are calculated from the transport and thermodynamic properties of the individual species 
provided in a recent study. The air model considered herein consists of 11 species: O 2 , N 2 , O, N, 
NO, 0+ N + , NO + , 0 ++ , N ++ , and e“. For pressures less than 10 -2 atm and temperatures 
greater than 14 000 K, the concentrations of N"^ 4- and O"^ become dominant. The differences 
in the computed values given herein and those obtained by other researchers may be attributed 
to the differences in the properties of the individual species used (including those due to the 
electron number density correction for the ionized-species transport properties), mixing laws 
used for the mixture, and to a lesser extent, the number of chemical species considered. 

The computed values of the various properties are curve fit as a function of temperature 
at constant values of pressure through the use of the least-squares curve-fit technique. A 
logarithmic interpolation relation is provided to obtain values at intermediate pressures. Further, 
a linear averaging procedure is suggested for the various curve-fit coefficients in order to obtain 
continuous values of the properties across the curve-fit boundaries. The curve fits are accurate 
to within 5 percent of the computed results for the pressure values and temperature range 
considered. These curve fits provide an efficient and fairly accurate means of computing the 
flow-field properties for equilibrium air calculations. 


NASA Langley Research Center 
Hampton, VA 23665-5225 
July 10, 1991 
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Appendix A 

Energy Equation in Terms of Frozen and Total Values of Specific Heat, 
Thermal Conductivity, and Prandtl Number 

The viscous-shock- layer form 3 of the energy equation for a chemically reacting multicom- 
ponent gas mixture can be written in terms of the frozen values of thermal conductivity and 
specific heat at constant pressure as follows (ref. 4): 

f u BT , BT\ ( u dp Bp\ 

P p -f\l + nK<9s V Bn) \1 + nnBs V Bn) 

= d_ ( K dT\ ( k , k &L 

dn \ f On ) \ 1 + uk r + n cos 6 ) * dn 


NS BT 

i= 1 


+ /x 



KU 

1 + UK 



(Al) 


If the species continuity equation (5) of reference 4 is multiplied by hj, summed over all the 
species, and then subtracted from equation (Al), the following form of the energy equation is 
obtained: 


p c Pf ( j 


_u BT BT\ _ f u. Bp 

+ nn Bs Bn J \ 1 + uk Bs 


, 9p\ 


l ~ 1 


K j COS 9 

+ [ + — 

1 + uk r + n cos 


du KU 


e) (*/£-£** 


+ M "ST ~ 


dn 1 + uk 


NS NS 

- > h iir L + v > h i~^ jL 

1 + UK 4-f ds dn 


2 — 1 


2 — 1 


(A2) 


Equation (A2) may be further rewritten, in terms of the total thermal conductivity and total 
specific heat at constant pressure, as follows: 


pCp 


u BT BT\ 

1- v — 

1 + uk Bs Bn J 


u Bp Bp 

1 + nn Bs Bn 


_B_ ( BT\ ( k 
Bn \ Bn ) \ 1 + 


j cos 0 
nn ' r + n cos 0 


+ 


K 


du 


KU \ 1 


+ ^ \dn 1 + nn) 


(A3) 


where Cp is the total specific heat at constant pressure and K is the total thermal conductivity 
defined in equations (1) and (2), respectively. The mass diffusion term in equations (Al) and 
(A2) is defined by equations (3), (4), or (5). 

3 The equations chosen for the purpose of illustration here are of the viscous-shock-layer form for the sake of simplicity. 
Similar treatment can be extended to the Navier-Stokes equations also. 
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Equations (A2) and (A3) do not contain the species production rate term ivj appearing 
in equation (Al). This is the term that is responsible for causing stiffness and convergence 
problems in the solution of equation (Al) when the flow is under (or near) chemical-equilibrium 
conditions. For such cases, forms of the energy equation given by equations (A2) and (A3) are 
more desirable. Unlike equation (A2), diffusion terms do not appear explicitly in equation (A3) 
but instead are contained in the definitions of total specific heat and total thermal conductivity 
(eqs. (1) and (2)). 

Reference 4 provides an alternate form of the energy equation (Al) in terms of total enthalpy 
for chemical-equilibrium flow calculations: 


(it 


OH dH_ 
nK ds dn 


dp KU^V 

~ v ^ l + pt~l — 

on 1 + nn 


dn 


+ 


(if 


j cos 9 

-I 

uk r + n cos 9 




(A4) 


where 


<f> = 


H dH 


P 


N PrJ 9n N Pr j ^ 


NS ac . 


dn 


NS 

P 'y ( h'i Ji T 

i—l 


P 




Pr ,/ 


("Pr ,/ - l) 


du 


u- 


{.IKll 


dn 1 + nn 


(A5) 


Similar to equations (A2) and (A3), equation (A4) also does not contain the species production 
rate term thj and is a desired form for equilibrium flow calculations. 

If the frozen Prandtl number Np r j appearing in equation (A5) is replaced by the total 
Prandtl number iVp r , equation (A5) becomes 


4 > = 


/i dh 
Np r dn 


+ pu 



KU 

1 + UK 


(A6) 


where the total Prandtl number Ap r is defined as (eq. (9)) 


Alpr — C p p/K 

and total values for C p and K are defined in equations (1) and (2), respectively. Equation (A6) 
is preferred over equation (A5) because the diffusion terms do not appear explicitly. 

Equations (A3) and (A4) (with $ defined in equation (A6)) have the same form as that for 
a perfect gas and may be solved in a similar manner without the need to solve for the species 
concentrations, provided the curve fits to the total values of Cp, K , and A [ > r are available for 
the ranges of temperatures and pressures of interest. 

There is a specific advantage to using equation (A3) instead of equation (A4) (with 4? given 
by equation (A6)) for equilibrium flow calculations because the expensive iterative process of 
obtaining the temperature (needed for computing the thermodynamic and transport properties) 
from the mixture enthalpy at every point in the flow field is avoided. Another advantage of using 
the energy equation with temperature as the dependent variable is that this form can also be 
used for finite-rate chemistry calculations. For such calculations, however, equation (A2), which 
separates the frozen and reactive components of C p and K, is more desirable. 
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Appendix B 

Sample Program To Evaluate Transport and Thermodynamic Properties for 
Equilibrium Air From Polynomial Curve Fit at a Given Pressure 

This appendix is a sample FORTRAN subroutine that evaluates the thermodynamic proper- 
ties for equilibrium air (with 1 1 constituent species) with polynomial curve fits used as functions 
of temperature. The same subroutine can also be modified for evaluation of the transport prop- 
erties. Different temperature ranges are used at different pressures. Properties evaluated near 
the temperature range boundaries are smoothed by linearly averaging the polynomial coefficients 
to assure continuous derivatives. The subroutine may be easily modified for different needs. 

SUBROUTINE THERMO (T,H) 

C Computes enthalpy of equilibrium air by approximating polynomials. 

C Polynomial coefficients are stored in arrays A1 to A5 and axe 
C linearly averaged at the temperature range boundaries. The 
C temperature ranges are for pressure at 0.0001 atm. 

C 

C input: T temperature, K 

C output: H enthalpy, kcal/g 

C 

DIMENSION A1 (6) ,A2(6) ,A3(6) ,A4(6) ,A5(6) 

DIMENSION P(5) ,C0EF(6,5) 

EQUIVALENCE (A1,C0EF) 

C 

C Coefficients are input for six temperature ranges 
C 


K=5 


L=6 

IF(T.GT. 11000.) GO TO 10 
K=4 


L=5 

IF(T. GT. 7000 . ) GO TO 20 
K=3 


L=4 

IF (T. GT.4500 . ) GO TO 30 
K=2 


L=3 

IF(T.GT. 2500. ) GO TO 40 
K=1 
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L=2 

PA=1 . 0 
PB=0 . 0 

IF(T. LE . 2000 . ) GO TG 50 
PB=0 . 002* (T-2000 . ) 

PA=1 . O-PB 
GO TO 50 
10 CONTINUE 
PA=0 . 0 
PB=1 . 0 

IFCT.GE. 18000.) GO TO 50 
PA=1 . 0 
PB=0 . 0 

IF (T.LE. 17500.) GO TO 50 
PB=0 . 002* (T-17500 . ) 

PA=1 . O-PB 
GO TO 50 
20 CONTINUE 
PA=1 . 0 
PB=0 . 0 

IFCT.LE. 10500.) GO TO 50 
PB=0 . 002* (T-10500 . ) 

PA=1 . O-PB 
GO TO 50 
30 CONTINUE 
PA=1 . 0 
PB=0 . 0 

IFCT.LE. 6500.) GO TO 50 
PB=0 . 002* (T-6500 . ) 

PA=1 . O-PB 
GO TO 50 
40 CONTINUE 
PA=1 . 0 
PB=0 . 0 

IFCT.LE. 4000.) GO TO 50 
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PB=0 . 002* (T-4000 . ) 

PA=1 . 0-PB 
50 CONTINUE 
C 

DO 100 1=1,5 

100 P (I)=PA*COEF (K , I)+PB*COEF (L , I) 

C 

TLN=ALOG(T/ 10000.) 

TLN2=TLN*TLN 

TLN3=TLN2*TLN 

TLN4=TLN3*TLN 

C 

H=P(1) *TLN4+P(2) *TLN3+P (3) *TLN2+P (4) *TLN+P (5) 
H=EXP (H) 

C 

RETURN 

END 
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Table L Curve-Fit Coefficients for Enthalpy of Equilibrium Air 
[Enthalpy obtained in kcal/g; curve fits obtained from enthalpies of individual species referenced to T = 0 K] 
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For temperatures less than 500 K, enthalpy for air may be obtained from the relation h = 0.24 x 10 


Table II. Curve-Fit Coefficients for Specific Heat of Equilibrium Air 
[Specific heat obtained in cal/g-K] 
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For temperatures less than 500 K, a value of 0.24 cal/g-K for the specific heat of air may be used. 



Table II. Concluded 
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For temperatures less than 500 K, a value of 0.24 cal/g-K for the specific heat of air may be used. 



Table III. Curve-Fit Coefficients for Compressibility of Equilibrium Air 
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For temperatures less than 500 K, a compressibility factor of unity may be used for air. 


Table IV. Curve-Fit Coefficients for Viscosity of Equilibrium Air 
[Viscosity obtained in poise] 
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Table V. Curve-Fit Coefficients for Thermal Conductivity of Equilibrium Air 
[Thermal conductivity obtained in cal/cm-sec-K] 
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For temperatures less than 500 K, Sutherland’s thermal conductivity law may be used: 
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Table VI. Curve-Fit Coefficients for Prandtl Number of Equilibrium Air 


a> 


i O O 

o 

o 

o 

o 

o 

o 

o 

O 

o 

o 

o 

o 

o 

o 

o o 

o 

o 

o 

o 

o 

o 

o o 

o 

o 

o 

o 

Q 

b 


ic 

to 

to 

to 

to 

to 

to 

o 

to 

to- 

to 

to 

to 

to 

o 

o 

to to 

to 

to 

to 

to 

o 

o 

to to 

to 

to 

to 

o 


P 

CM 

b 

b 

CN 

b 

b 

CN 

o 

CN 

b- 

CN 

CN 

b- 

CN 

to 

o 

lb CN 

CN 

t- 

CN 

CN 

to 

o 

b CN 

b 

b 

CN 

tO 

o 



CM 

CO 

to 

GO 

o 

Tp 

oo 

to 

CN 



o 

CN 

b- 

o 

00 

CN tO 

00 

1— 1 

TP 

00 

CO 

o 

CN tO 

b 

CO 

00 

to 

o 

e 





rH 

i— i 

H 

CN 




iH 

i— l 

rH 

CN 

CN 



1—1 

i— l 

t-H 

CN 

co 



t-H 

r-H 

CN 

CO 

<D 

ft 

a 

o 

bC 





P 

o 

to 

to 

tO 

to 

tO 

to 

to 

O 

tO 

to 

to 

tO 

tO 

to 

o 

O to 

LO 

to 

to 

to 

to 

o 

O to 

LO 

to 

to 

to 

o 

d 

to 

CN 

b 

b 

CN 

b 

b 

CN 

to 

CN 

b- 

CN 

CN 

i- 

CN 

to 

tO b* 

CN 

CN 

lb 

CN 

CN 

to 

to b 

CN 

b 

b 

CN 

to 

r cu 

j-t 


CN 

CO 

to 

00 

o 

Tp 

00 


CN 

H 


O 

CN 

b- 

o 

CN 

to 

00 

rH 

Tp 

00 

CO 

CN 

to 

b 

CO 

GO 

to 

H 







1— 1 

iH 

rH 





rH 

rH 

H 

CN 




i— ! 

rH 

rH 

CN 




t-H 

r-H 

CN 



H 

o 

o 

CN 

rH 

CO 

CO 

CD 

H 

rH 

iH 

H 

CN 

CO 

CO 

to 

rH rH 

rH 

CN 

CN 

Tp 

tP 

to 

CN i~H 

i-H 

tP 

TP 

tP 

CO 



o 

t 

o 

+ 

q 

o 

4- 

q 

O 

1 

o 

[ 

o 

1 

o 

1 

O 

i 

o 

1 

o 

o 

1 

o 

1 

O 

1 

o 

I 

o 

I 

o 

1 

o o 

i i 

o 

i 

O 

j 

O 

1 

O 

I 

O 

j 

o 

1 

o o 

I 1 

o 

1 

O 

1 

o 

1 

O 

1 

o 

J 



I 

63 

1 

63 

1 

63 

1 

63 

1 

63 

1 

63 

1 

63 

1 

63 

1 

63 

1 

63 

I 

63 

1 

63 

1 

63 

1 

63 

i 1 

63 63 

i 

63 

1 

63 

1 

63 

1 

63 

1 

63 

1 

63 

I 1 

63 63 

1 

63 

1 

63 

1 

63 

1 

63 

1 

63 



CO 

H 

o 

co 

b 

00 

00 

to 

CN 

b- 

b~ 


CO 

CO 

rH 

to 

05 O 

b- 

00 

o 

05 

CN 

CN 

CO co 

b 

o 

LO 

o 

CN 


n i 

Tp 

CO 

CN 

CO 

o 

CN 

CN 

o 

CD 

CD 

CD 

r- 

CO 

s 

05 

to 

CO O 


CN 

o 

co 

CO 

Tp 

CN CO 

CD 

to 

05 

b 

05 



00 

T— 1 

05 

CO 

b 

to 

rH 

CN 

tO 

H 

O 

tr- 

rH 

o 

CD 

05 

O CN 

00 

rH 

CD 

to 

CD 

05 

O 05 

oo 

00 

b 

o 

TP 



Tp 

o 

b 

to 

Tf 

rH 

rH 

o 

H 

CO 

b- 

H 

H 

o 

co 

CD 

to 

rH 

co 

i“H 

to 

05 

CO 

rH CD 

to 

CD 

rH 

CD 

CO 



oo 

co 

Tp 

CN 

b 

05 

00 

Tp 

CO 

tO 

oo 

CD 

CN 

to 

05 

rH 

H 05 

s 

05 

CN 

CD 

3 

b 

tP CN 

oo 

co 

co 

t-H 

00 



CN 

Tp 

b 

CD 

rH 

CN 

00 

b 

b- 


b- 

CN 

CD 

rH 

rH 

CN 

tO CO 

lb 

rH 

CO 

CN 

tO 05 

o 

00 

CO 

b 

CN 



05 

rH 

CN 

to 

00 

CN 

Tp 

co 

to 

05 

05 

CN 

CO 

rH 

05 

CO 

CN ^ 

CO 

i-H 

co 

CO 

Tl 4 

CD 

CD to 

CO 

r-H 

CN 

CN 

co 



d 

o 

1 

d 

d 

1 

o 

d 

1 

o 

d 

d 

d 

1 

d 

1 

d 

d 

o 

1 

d 

1 

d 

d d 

I 

o 

1 

d 

o 

o 

1 

o 

O 

d d 

I 

o 

d 

d 

d 

d 



o 

rH 

i-H 

o 

rH 

rH 

H 

to 

o 

rH 

rH 

Q 

o 

CN 

rH 

co 

o o 

rH 

rH 

Q 

CN 

CN 

CO 

rH rH 

o 

CN 

CO 

CN 

rH 



o 

O 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

O 

o 

o 

o o 

o 

o 

o 

o 

O 

o 

o o 

o 

O 

o 

o 

o 



4- 

+ 

+ 

+ 

+ 

1 

1 

1 

+ 

+ 

+ 

+ 

+ 

1 

1 

i 

+ + 

4- 

1 

4* 

1 

1 

1 

i + 

+ 

1 

1 

i 

1 



63 

63 

63 

63 

63 

q 

q 

63 

63 

63 

63 

63 

63 

63 

63 

63 

63 63 

q 

63 

q 

63 

63 

63 

63 63 

63 

63 

q 

q 

63 



to 

b 

CO 

CD 

tP 

to 

CO 

CN 

o 

oo 

to 

05 


rH 

b- 

CN 

CO CO 

o 

05 

CN 

rH 

b 

CO 

00 tP 

CN 

co 

CN 

Tp 

Tp 


ft 

b 

CD 

b 

b 

tP 

CN 

b- 

to 

CD 

o 

b- 


CO 

b- 

to 

CN 

CO CO 

o 

CD 

CN 

CD 

CO 

rH 

i—i co 

to 

CN 

Tp 

b 

CD 


63 

CO 

b 

to 

to 

CN 

to 

05 

b 

tO 

CD 

05 

00 

o 

H 

o 

to 

05 05 

to 

CD 

rH 

O 

Tf 4 

rH 

b o 

b 

O 

05 

CN 

tP 


CN 

O 

CO 

rH 

05 

to 

rH 

CD 

b- 

CN 

CO 

00 

o 

rH 

to 

to 

CN 

rH 

CO 

o 

00 

CD 

o 

CD CD 

b 

CO 

TP 

CO 

05 



b 

to 

CO 

o 

CD 

CD 

o 

CO 


rH 



co 

oo 

05 

oo 

05 tO 

CN 

T* 

o 

oo 

rH 

CD 

CN tT 

b 

CN 

00 

b 

05 



CN 

CO 

CO 

CO 

O 

tO 

b 

b- 

00 

b- 

05 

rH 

05 

b- 

b- 

CD 

CN b- 

CN 

to 

o 


o 

CO 

CO CN 

CD 

TP 

CN 

05 

Tp 



Tf 

CN 

CD 

rH 

Tf 

rH 

co 

CD 

CN 

rH 

CN 

05 

iH 

00 

00 

CO 

H 00 

rH 

CD 

CN 

to 

tP 

00 

r-H i-H 

05 

tO 

to 

CN 

Tp 



d 

i 

o 

o 

1 

o 

d 

1 

o 

d 

1 

O 

o 

1 

o 

O 

d 

1 

d 

1 

d 

o 

d 

i 

d d 

1 

d 

O 

I 

o 

I 

o 

o 

I 

d 

1 

d d 

I 

d 

1 

o 

d 

1 

d 

1 

d 

1 



© 

CN 

CN 

rH 

CN 

o 

rH 

CN 

o 

CN 

CN 

CN 

rH 

o 

rH 

rH 

O i — < 

CN 

o 

rH 

o 

o 

rH 

O CN 

CN 

o 


o 

rH 



o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

O 

o 

o 

o 

o 

o 

o o 

o 

o 

o 

o 

o 

o 

o o 

O 

o 

o 

o 

o 



+ 

+ 

+ 

+ 

+ 

+ 

+ 

1 

+ 

+ 

+ 

+ 

+ 

+ 

4- 

1 

4* 4- 

+ 

+ 

+ 

+ 

+ 

1 

+ + 

+ 

+ 

1 

4- 

+ 



63 

63 

63 

63 

63 

63 

63 

63 

63 

63 

63 

63 

63 

63 

63 

63 

63 63 

63 

63 

63 

63 

63 

63 

6)63 

63 

63 

63 

63 

63 



i-H 

to 

Tf 

Tp 

b 

CN 

b 

b 

co 

b- 

00 


CN 

CN 

o 

rH 

b- b- 



b 

co 

CO 

CN 

b CD 

b 

to 

05 

05 

CO 


ft 

to 

CN 

CO 

rH 

rH 

tO 

tP 

b- 

rH 

CN 

o 

05 

tO 

o 

rH 

CD 

b- 00 

CN 

rH 

CN 

CN 

CD 

tP 

00 o 

b 

to 

O 

tO 

o 


C| 

b 

<75 

Tp 

b 

tO 


CO 

CO 

rH 

to 

H 

00 

to 

b- 



05 O 

rH 

CD 

Tt 4 

CN 

CO 

CN 

00 CO 

b 

to 

05 

CO 

b 


b 

O 

on 

rH 

05 

tH 

05 

CN 

CO 

H 4 

CN 

o 

o 

b- 

oo 

rH 

o to 

CO 

to 

CN 

to 

rH 

O 

05 to 

00 

o 

00 

o 

co 



CO 

o 

CD 

05 

CD 

tO 

O 

oo 


O 

O 

CD 

o 

CD 

to 

LO 

00 to 

CN 

b- 

00 

00 

tO 

CO 

1— i CN 

05 

05 

rH 

05 

CD 



CN 

to 

b 

H 

O 

CN 

rH 

b 

00 

CN 

to 

H 

o 

CD 

co 

CD 

^ b~ 

CD 

CO 

o 

00 

Tf 

Tp 

O O 

rH 

to 

Tp 

CN 

Tp 



to 

rH 

to 

rH 

OO 

Tp 

rH 

CN 

CO 

rH 

CO 

H 


CN 

co 

rH 

rH cD 

rH 

00 

to 

i-H 

rH 

TP 

rH rH 

r-H 

CN 

CO 

rH 

CN 



d 

o 

i 

o 

o 

1 

© 

d 

f 

d 

d 

1 

d 

d 

1 

o 

1 

d 

d 

o 

i 

d 

1 

o 

d d 
1 

o 

j 

o 

o 

o 

1 

o 

O 

o d 
1 1 

d 

d 

1 

d 

1 

d 

d 



CN 

CN 

CO 

rH 

CO 

rH 

CN 

o 

rH 

CN 

CO 

co 

CN 

rH 

CN 

o 

H CN 

co 

rH 

CN 

rH 

rH 

rH 

O CN 

CN 

t-H 

rH 

rH 

CN 



O 

O 

o 

o 

o 

o 

o 

o 

o 

O 

o 

o 

o 

o 

O 

o 

o o 

o 

o 

o 

o 

o 

o 

o o 

o 

o 

o 

o 

O 



1 

+ 

4 

+ 

4 

+ 

+ 

+ 

1 

+ 

+ 

+ 

+ 

+ 

-4 

+ 

1 4- 

4 

+ 

+ 

+ 

+ 

+ 

+ + 

4- 

4 

+ 

+ 

+ 



q 

63 

63 

63 63 

q 

q 

q 

^q 

63 

cq 

cq 

sq 

cq 

cq 

cq 

qq 

63 

63 

63 

q 

q 

q 

qq 

q 

63 

q 

q 

q 



o 

o 

rH 

CN 

to 

05 

o 

05 

CN 

<o 

CN 

o 

co 

ot 4 

GO 

co 

to 00 

C5 

rH 

rH 

CD 

CN 

CO 

O CN 

CN 

o 

CN 

CN 

b 


ft 

o 

rH 

00 

rH 

rH 

05 

b 

CO 

H 

rH 

CN 

CN 

CN 

o 

b- 

to 

oo b- 

b~ 

to 

to 

05 

to 

oo 

oo to 

CN 

r-H 

O 

b 

00 



CN 

to 

to 

CO 

05 

tO 

CD 

05 

tO 

00 

CD 

NT 

oo 

H 


o 

CO 1— ( 

00 

r— 1 

tP 

CN 

05 

CD 

Tp TP 

to 

i— 1 

CD 

co 

00 


H 

CN 

Tf 

rH 

CO 

CN 

05 

Tp 

00 

rH 

05 

CO 


to 

b- 

CD 

H* 05 

^f 4 

05 

00 

05 

TP 

to 

00 b 

o 

b 

05 

to 

b 



b 

TP 

o> 

co 

CD 

tP 

CO 

to 


o 

tO 

to 

00 

CO 

05 

00 

rH rH 

CD 

00 

00 

CN 

to 

CD 

CD i“H 

CO 

CN 

rH 

b 

CO 



CO 

CD 

LO 

o 

05 

b 

CD 

CN 

tO 

H 

O 

rH 

05 

o 

CO 

tO 

o to 

o 

o 

co 

CN 

to 

rH 

i— 1 i-H 

CO 

CN 

tO 

b 

b 



00 

tP 

CN 

CD 

b 

to 

rH 

H 

CO 


CN 

rH 

CO 

oc 

CD 

CO 

00 CN 

1— 1 


CD 

CO 

CN 

rH 

CO tP 

b 

TP 

rH 

CN 

CD 



d 

d 

d 

1 

d 

d 

1 

d 

d 

1 

d 

d 

1 

d 

d 

o 

1 

d 

! 

d 

d 

o 

I 

d d 

d 

d 

i 

d 

1 

d 

d 

1 

d 

1 

d d 

d 

i 

d 

d 

d 

! 

d 

1 



o 

CN 

co 

CN 

tP 

CN 

co 

rH 

o 

CN 

co 

co 

CO 

CN 

CO 

rH 

O CN 

co 

rH 

CO 

CN 

CN 

CN 

O CN 

CO 

CN 

CN 

CN 

co 



o 

O 

o 

o 

O 

O 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o o 

o 

o 

o 

O 

O 

O 

o o 

o 

O 

O 

o 

o 



+ 

+ 

4- 

+ 

+ 

+ 

+ 

+ 

+ 

+ 

+ 

+ 

+ 

4- 

+ 

+ 

+ + 

4- 

+ 

4- 

+ 

+ 

+ 

4- 4- 

+ 

4- 

4- 

+ 

+ 



q 

q 

q 

63 63 

q 

q 

£q 

eq 

eq 

cq 

tq 

tq 

q 

tq 

q 

qq 

q 

q 

q 

q 

q 

q 

qq 

q 

q 

q 

q 

q 



CD 

CN 

i—i 

CO 

CD 

b- 

to 

CN 

o 

b- 

b- 

oo 

CN 

o 

05 

b- 

CD tO 

00 

ib 

TP 

CO 

CD 

b 

CD 05 

o 

o 

00 

co 

b 


ft 

tO 

00 

tp 

05 

Tp 

05 

Tp 

00 

to 


CD 

00 

00 

00 

tO 

CN 

00 to 



o 

r-H 

05 

TP 

CD 00 

Tp 

CN 

CD 

b 

00 



i— l 

Tp 

tp 

to 

CD 

CN 

Tp 

H 

o 

rH 

co 

b- 

uo 

05 

H 4 

Nt 4 

CD 00 

o 

lb 

o 

i-H 

rH 

TP 

tO 00 

CD 

CN 

00 

to 

LO 


b 

to 

to 

CO 

00 

H 

CO 

05 

05 

rH 

H 

tO 

CO 

H 

00 

tO 

05 H 

b- 

CO 

o 

i-H 

oo 

o 

i“H CN 

Tp 

05 

oo 

05 

CO 



00 

Tf 

CD 

CN 

CO 

CD 

05 

o 

to 

to 

00 

H 

CD 

H 

b- 

rH 

CN CO 

3 

CN 

CD 

CD 

b 

CN 

i — i lO 

t-H 

CN 

CO 

to 

o 



TP 

05 

b 

05 

05 

00 

rH 

CN 

CD 

b- 

05 

to 

00 

O 

o 

b- 

CD tO 

CO 

05 

b 

rH 

to 

CD O 

CN 

05 

rH 

05 

CN 



CN 

CD 

tO 

1— t 

CO 

CO 

rH 

CN 

rH 

CD 

to 

rr 

rH 

CO 

CD 

co 

rH T? 

CO 

CN 

co 

CN 

CN 

rH 

CN 00 

CN 

CN 

CN 

CN 

05 



d 

1 

d 

1 

d 

d 

1 

d 

d 

_L 

d 

d 

1 

d 

! 

d 

J_ 

d 

1 

d 

d 

d 

_L 

d 

f 

d 

d d 
1 1 

d 

d 

d 

d 

1 

d 

d 

d d 
1 1 

d 

d 

1 

d 

! 

d 

o 



o 

CN 

co 

CN 

H 

CO 

co 

IM 

o 

CN 

CO 

co 

co 

CN 


CN 

O CN 

CO 

CN 

CO 

CN 

CN 

CN 

O CN 

CO 

CN 

CO 

CO 

TP 



o 

O 

o 

o 

o 

o 

o 

O 

o 

O 

o 

o 

o 

O 

o 

o 

O O 

o 

o 

o 

O 

O 

O 

O O 

o 

o 

o 

o 

O 



+ 

-F 

+ 

4- 4- 

+ 

+ 

+ 

+ 

+ 

+ 

+ 

+ 

4- 

4- 

+ 

4- 4- 

+ 

4- 

4- 

+ 

4- 

+ 

4- + 

4- 

+ 

+ 

4- 

+ 



q 

q 

q 

63 63 

q 

q 

63 

cq 

cq 

cq 

tq 

tq 

tq 

q 

q 

qq 

q 

q 

q 

q 

q 

q 

qq 

q 

q 

q 

q 

q 


. 

00 

CD 

b 

Tp 

CN 

to 

o 

tH 

H 

00 

oo 

05 

00 

b- 

to 

CO 

CN tO 

to 

05 

05 

b 

S 

00 

b b 

00 

CD 

rH 

oo 

to 


ft 

1—1 

CN 

o 

H 

CN 

CN 

CO 


b- 

o 

b- 

00 

CN 

CN 


rH 

rH CN 

o 


oo 

rH 

rH 

b 00 

CD 

05 

b 

CD 

b 



CO 

05 

05 

b 

CD 

00 

to 

GO 

o 

CD 

to 

05 

co 

CN 

CN 


o oo 

rH 

o 

b 

b 

to 

b 

CO co 

00 

05 

b 

rH 

05 


to 

H 

o 

CN 

Tf 

Tp 

Tp 

GO 

CO 

to 

b- 

H 

b- 

CO 



CD CD 

05 


05 

CD 

CO 

CD 

00 o 

i— i 

rH 

Tf 

CN 

CN 



o> 

oo 

CD 

b 

CN 

co 

rH 

05 

oo 

to 

b- 

CN 

o 

to 

O 

CN 

o o 

o 


co 

CO 

CN 

Tp 

Tp 05 

tJ 4 

to 

CO 

to 

CN 



CD 

o 

o 

rH 

CD 

o 

Tf 

co 


CO 

00 

O 

H 

rH 

CO 

H 4 

^ CN 

TJ 4 


CD 

tP 

r-H 

00 

lO rH 

CD 

Tp 

o 

CN 

O 



b 

Tf 

to 

co 

b 

iH 

CO 

H 

b- 

rf 

CD 

b- 

co 

05 

CN 

iH 

b- CO 

Tt 4 

CN 

05 

05 

b 

b 

b CD 

CN 

b 

rH 

i-H 

tO 



d 

o 

d 

X 

o 

d 

1 

o 

d 

1 

o 

d 

o 

d 

d 

1 

d 

1 

d 

d 

d 

1 

d d 

d 

d 

d 

1 

d 

d 

d 

d d 

d 

1 

d 

d 

d 

1 

d 

1 

of 

(H 

P 

P 

N* 

| 








C0 








<N 







H 






c n 

P 

o 








1 

O 








1 

o 







1 

o 






o 

u 

0^ 

d 

T— ( 








tH 








rH 







rH 







30 


a For temperatures less than 500 K, the following relation (employing the Sutherland values for the viscosity and thermal conductivity) may 
be used to obtain Prandtl number for air: 

Vp r = 0.24 fJ, s /K s 
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Figure 2. Enthalpy of equilibrium air. 
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Figure 2. Continued. 



Figure 2. Concluded, 




Figure 3. Specific heat at constant pressure of equilibrium air. 



42 


Figure 3. Continued. 
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Figure 3. Concluded. 



Figure 4. Compressibility factor for equilibrium ai 



Figure 4. Continued. 
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Figure 4. Concluded. 
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Figure 5. Continued. 
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Figure 5. Concluded. 
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Figure 6. Thermal conductivity for equilibrium ai 



Figure 6. Continued. 



Figure 6. Concluded. 



Figure 7. Prandtl number for equilibrium air. 
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Figure 7. Continued. 
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Figure 7. Concluded. 
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(b) Total specific heat at constant pressure. 
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(f) Total Prandtl number. 
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Figure 9. Comparison of dimensionless enthalpy values for equilibrium air. 




Figure 10. Concluded. 








Figure 12. Comparison of compressibility factor values from reference 7 and present computation for equilib- 
rium air. 




Figure 13. Comparison of viscosity values from several sources for equilibrium air. 
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Figure 14, Concluded, 






Figure 15. Concluded. 







